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Abstract — Recently, quantum error-correcting codes were pro- 
posed that capitalize on the fact that many physical error models 
lead to a significant asymmetry between the probabilities for 
bit flip and phase flip errors. An example for a channel which 
exhibits such asymmetry is the combined amplitude damping and 
dephasing channel, where the probabilities of bit flips and phase 
flips can be related to relaxation and dephasing time, respectively. 
We give systematic constructions of asymmetric quantum stabi- 
lizer codes that exploit this asymmetry. Our approach is based 
on a CSS construction that combines BCH and flnite geometry 
LDPC codes. 

L Introduction 

In many quantum mechanical systems the mechanisms for 
the occurrence of bit flip and phase flip errors are quite 
different. In a recent paper loffe and Mezard [10] postulated 
that quantum error-correction should take into account this 
asymmetry. The main argument given in [10] is that most of 
the known quantum computing devices have relaxation times 
(Ti) that are around 1 — 2 orders of magnitude larger than 
the corresponding dephasing times (T2). In general, relaxation 
leads to both bit flip and phase flip errors, whereas dephasing 
only leads to phase flip errors. This large asymmetry between 
Ti and T2 suggests that bit flip errors occur less frequently 
than phase flip errors and a well designed quantum code 
would exploit this asymmetry of errors to provide better 
performance. In fact, this observation and its consequences for 
quantum error correction, especially quantum fault tolerance, 
have prompted investigations from various other researchers 
[1], [8], [20]. 

Our goal will be as in [10] to construct asymmetric quantum 
codes for quantum memories and at present we do not consider 
the issue of fault tolerance. We first quantitatively justify how 
noise processes, characterized in terms of Ti and T2, lead 
to an asymmetry in the bit flip and phase flip errors. As a 
concrete illustration of this we consider the amplitude damping 
and dephasing channel. For this channel we can compute the 
probabilities of bit flip and phase flips in closed form. In 
particular, by giving explicit expressions for the ratio of these 
probabilities in terms of the ratio Ti /T2, we show how the 
channel asymmetry arises. 

After providing the necessary background, we give two 
systematic constructions of asymmetric quantum codes based 
on BCH and LDPC codes, as an alternative to the randomized 
construction of [10]. 



II. Background 
Recall that a quantum channel that maps a state p to 

(1 -Px-Py -Pz)p+PxXpX-Fp^jYpY-Fp^ZpZ, (1) 

withl= [JO], X=[oi],Y=[o-^],Z=[J_°] is called 
a Pauli channel. For a Pauli channel, one can respectively 
determine the probabilities Px,Pv>Pz that an input qubit in 
state p is subjected to a Pauli X, Y, or Z error 

A combined amplitude damping and dephasing channel £ 
with relaxation time Ti and dephasing time T2 that acts on a 
qubit with density matrix p ~ (Pij )i,je{o,i } for a time t yields 
the density matrix 
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This channel is interesting as it models common decoherence 
processes fairly well. We would like to determine the proba- 
bility Px, Py, and pz such that an X, Y, or Z error occurs 
in a combined amplitude damping and dephasing channel. 
However, it turns out that this question is not well-posed, since 
£ is not a Pauli channel, that is, it cannot be written in the 
form ([T]i. However, we can obtain a Pauli channel £j by a 
technique called twirling [7], [5]. In our case, the twirling 
consists of conjugating the channel £ by Pauli matrices and 
averaging over the results. The resulting channel £j is called 
the Pauli-twirl of £ and is explicitly given by 

1 



£.[9) = 4 



Y_ At£(ApAt)A. 
Ae{i,x,Y,z; 

Theorem 1: Given a combined amplitude damping and 
dephasing channel £ as above, the associated Pauli-twirled 
channel is of the form 



[1 



c -Py -Pz)P + PxXpX + pyYpY + PzZpZ, 

where p^ = p^ = (1 - e"*/^i )/4 and p^ = 1/2 - px - 
jC^*/^^. In particular, 

n 1 - pt/T, (1-T,/T2) 

^ = 1+2- 



Px eVT, _ 1 

If t <C Ti , then we can approximate this ratio as 2Ti /T2 — 1 . 
Proof: The Kraus operator decomposition [18] of £ is 
2 

£:(p) =^AkpA[, (2) 

k=0 



where An 



1 

\/l-A-Y 
t/Ti 



;Ai 



and VI — Y — A — e 
the Kraus operators At as 

1 



Y 





VX 

p-t/T 



•A, = 
' . We can rewrite 



Ao = 



Ai 



VA, V\ 



Z, A2 



-X 



2 

2 ' 2 2 2i 

Rewriting £(p) in terms of PauU matrices yields 

■2\/T 
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^IpZ - ^Zpl + ^XpY - ^YpX. 
4 4 4i 4i 



(3) 



It follows that the Pauli-twirl channel £j is of the claimed 
form, see [5, Lemma 2]. Computing the ratio Pz/Px we get 



Vx 
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et/T, _ I 

If t ^ Ti , then we can approximate the ratio as 2Ti /T2 — 1 , 
as claimed. ■ 

Thus, an asymmetry in the Ti and T2 times does translate 
to an asymmetry in the occurrence of bit flip and phase flip 
errors. Note that Px = Py indicating that the Y errors are as 
unlikely as the X errors. We shall refer to the ratio Pz/px as 
the channel asymmetry and denote this parameter by A. 

Asymmetric quantum codes use the fact that the phase flip 
errors are much more likely than the bit flip errors or the 
combined bit-phase flip errors. Therefore the code has different 
error correcting capability for handling different type of errors. 
We require the code to correct many phase flip errors but it is 
not required to handle the same number of bit flip errors. If we 
assume a CSS code [4], then we can meaningfully speak of X- 
distance and Z-distance. A CSS stabilizer code that can detect 
all X errors up to weight dx — 1 is said to have an X-distance 
of dx- Similarly if it can detect all Z errors upto weight dz — 1, 
then it is said to have a Z-distance of dz. We shall denote such 
a code by [[n, k, dx/ dz]] q to indicate it is an asymmetric code, 
see also [19] who was the first to use a notation that allowed to 
distinguish between X- and Z-distances. We could also view 
this code as an [[n, k, min{dx, dz}]]q stabilizer code. Further 
extension of these metrics to an additive non-CSS code is an 
interesting problem, but we will not go into the details here. 

Recall that in the CSS construction a pair of codes are 
used, one for correcting the bit flip errors and the other for 
correcting the phase flip errors. Our choice of these codes 
will be such that the code for correcting the phase flip errors 
has a larger distance than the code for correcting the bit flip 
errors. We restate the CSS construction in a form convenient 
for asymmetric stabilizer codes. 



Lemma 2 (CSS Construction [4]): Let Cx,Cz be linear 
codes over with the parameters [n,kx]q, and [n,kz]q 
respectively. Let C Cz. Then there exists an [[n, kx + 
kz — n, dx/dz]]q asymmetric quantum code, where dx = 
wt(Cx\C^) and dz =wt(Cz\C;^). 

If in the above construction dx = wt(Cx) and dz = wt(Cz), 
then we say that the code is pure. 

In the theorem above and elsewhere in this paper Fq denotes 
a finite field with q elements. We also denote a q-ary narrow- 
sense primitive BCH code of length n = q"^ — 1 and design 
distance 6 as BCTi.{S). 

III. Asymmetric Quantum Codes from LDPC Codes 

In [10], loffe and Mezard used a combination of BCH 
and LDPC codes to construct asymmetric codes. The intu- 
ition being that the stronger LDPC code should be used for 
correcting the phase flip errors and the BCH code can be 
used for the infrequent bit flips. This essentially reduces to 
finding a good LDPC code such that the dual of the LDPC 
code is contained in the BCH code. They solve this problem 
by randomly choosing codewords in the BCH code which are 
of low weight (so that they can be used for the parity check 
matrix of the LDPC code). However, this method leaves open 
how good the resulting LDPC code is. For instance, the degree 
profiles of the resulting code are not regular and there is little 
control over the final degree profiles of the code. Furthermore, 
it is not apparent what ensemble or degree profiles one will 
use to analyze the code. 

We propose an alternate scheme that uses LDPC codes 
to construct asymmetric stabilizer codes. We propose two 
families of quantum codes based on LDPC codes. In the 
first case we use LDPC codes for both the X and Z channel 
while in the second construction we will use a combination of 
BCH and LDPC codes. But first, we will need the following 
facts about generalized Reed-MuUer codes and finite geometry 
LDPC codes. 

A. Finite Geometry LDPC Codes ([14], [21]) 

Let us denote by EG(m,p'') the Euchdean finite geometry 
over Fps consisting of p"^* points. For our purposes it suffices 
to use the fact that this geometry is equivalent to the vector 
space F^s . A |j.-dimensional subspace of F^s or its coset is 
called a ^-flat. Assume that < p-i < 1^2 < Ta. Then we 
denote by Neg(H2, Hi , s,p) the number of |a.i -flats in a 10.2 - 
flat and by AegIttl, \i2, M-i , s,p), the number of |a.2-flats that 
contain a given |j,i-flat. These are given by (see [21]) 



Neg(M-2,M-i,s,p) = q 



11 nM-i -1- 



i=l 



qra-i+l _ 1 

Aeg(ttl, M.2, M-1 ,s,p) = Yl q^2-i+i _ 1 ' 

1=1X1+1 ^ 



,(4) 



(5) 



where q = p**. Index all the |i.i -flats from i = 1 to n = 
MEG(Tn, M-i , s,p) as F^. Let F be a (a.2-flat in EG(m,p'^). Then 



we can associate an incidence vector to F with respect to the 
\X} flats as follows. 

J . I ij = 1 if Fj is contained in F 1 
1 ' ij = otherwise. J ' 

Index the |a,2-flats from j = 1 to J = Neg(tti, |j.2, s,p). 
Construct the J x n matrix Hg^ (m, |X2, M-i , s,p] whose rows 
are the incidence vectors of all the |j.2 -flats with respect to the 
P-i -flats. This matrix is also referred to as the incidence matrix. 
Then the type-I Euclidean geometry code from |a^2 -flats and 
(j-i -flats is defined to be the null space, i. e., Euclidean dual 
code) of the Fp-Hnear span of Hg^ (m, |j.2, M-i , s,p). This is 
denoted as (m, |j.2 , M-i , s , p ) . Let Hj,g' (m, |j.2 , M-i , s , p ) = 
^EG l-Li , s,p)*. Then the type-II Euclidean geometry 

code CgQ (m, IJ.2, H-i , s,p) is defined to be the null space of 
HgQ (m, |J.2, M^i , s,p). Let us now consider the p.2-flats and 
(j-i -flats that do not contain the origin of EG(m,p''). Now 
form the incidence matrix of the |j,2 -flats with respect to 
the p-i -flats not containing the origin. The null space of this 
incidence matrix gives us a quasi-cyclic code in general, which 
we denote by CgQ (m, |X2, M-i , s, p), see [21]. 

B. Generalized Reed-Muller codes ([12]) 

Let a be a primitive element in Fqm . The cyclic generalized 
Reed-Muller code of length q — 1 and order -v is defined as 
the cyclic code with the generator polynomial whose roots a' 
satisfy < j < ra(q — 1)— "v— 1. The generalized Reed-Muller 
code is the singly extended code of length q"^. It is denoted 
as GRMq(-i/,m). The dual of a GRM code is also a GRM 
code [2], [3], [12]. It is known that 

GRMq(-i/,m,)-L =GRMq(^-L,ra), (6) 

where -y^ = Ta(q — 1) — 1— "v. 

Let C be a linear code over F^s . Then we define C\r^, the 
subfield subcode of C over Fq as the codewords of C which 
are entirely in FJJ, (see [9, pages 116-120]). Formally this can 
be expressed as 

C\w, = {c e C I c e F^}. (7) 

Let C C ¥^i. The the trace code of C over Fq is defined as 

trqVq(C) ={trqvq(c) Ice C}. (8) 

There are interesting relations between the trace code and the 
subfield subcode. One of which is the following result which 
we will need later. 

Lemma 3: Let C C FJJi. Then CIf,, the subfield subcode 
of C is contained in trqi/q(C), the trace code of C. In other 
words 

C|f, C trqi/q{C). 

Proof: Let c G C|f, C F];^ and a G Fqi. Then 
trqi /q (ac) = c trqi/q (oc) as c e Fq . Since trace is a surjective 
form, there exists some a E Fqi, such that trqi/q(a) = 1. 
This implies that c £ trqi/q (C). Since c is an arbitrary element 
in C\w^ it foflows that C|f, C trqi/q(C). ■ 



Let q = p*, then the Euclidean geometry code of order r 
over EG(m.,p*) is defined as the dual of the subfield subcode 
of GRMq((q - l)(ra-r- 1),m), [3, page 448]. The type-I 
LDPC code C^q' (m, ix, 0, s,p) code is an Euclidean geometry 
code of order 1 over EG(Ta,p'*), see [21]. Hence its dual 
is the subfield subcode of GRMq((q — l)(Ta — M.),tu) code. 
In other words, 

c|,[;'(m,H,0,s,p)^ =GRMq((q-l)(m-^i),m)|F,. (9) 
Further, Delsarte's theorem [6] tells us that 
^gV.^.O.s.p) = GRMq((q-1)(m-H),m)|^^, 

- trq/p (GRMq((q-1)(m-H),m)^) 
= trq/p(GRMq(^(q-1)-1,m)). 

Hence, Cg^ (m, 0, s, p) code can also be related to 
GRMq(M.(q - 1) - l,ra) as 

4y(m,H,0,s,p) =trq/p(GRMq(H(q-1)-l),m). (10) 

C. New families of asymmetric quantum codes 

With the previous preparation we are now ready to construct 
asymmetric quantum codes from finite geometry LDPC codes. 

Theorem 4 (Asymmetric EG LDPC Codes): Let p be a 
prime, with q = p'' and s > 1,Ta > 2. Let 1 < \iz. < 
and m — + 1 < M^x < tti. Then there exists an 

[[p'"^lcx + k,-p'-^dx/dJ]p 

asymmetric EG LDPC code, where 

kx = dimCgc (m, Hx,0, s,p); k^ = dim (m, Hz, 0, s,p). 

For the distances dx > AEofTU, |a.x , M-x ^ 1 , s, p) + 1 and > 
AegIttl, |J.z, M-z - 1,s,p) + 1 hold. 

Proof: Let = C^q (m, |j.z,0, s,p). Then from equa- 
tion ([Tol l we have 

Cz = trq/p(GRMq(^z(q-1)-l,m). 

By Lemma |3] we know that 

C, D GRMq(^z(q-^)-^,m)|F,, 

C, D GRMq((q-1)(m-(m-Hz + 1)),m)|F,, 

where the last inclusion follows from the nesting property of 
the generalized Reed-Muller codes. For any order p-x such that 
m - lo^z + 1 < M-x < TTL, let Cx = CgQ (ra, |j.x, 0, s, p). Then 
Cx is an LDPC code whose dual =GRMq((q - l)(Ta- 
M-x),Tn)lFp is contained in C^. Thus we can use Lemma[2]to 
form an asymmetric code with the parameters 

[[p'"MCx+kz-p'-«,dx/dJ]p 

The distance of and Cx are at lower bounded as dx > 
AEG(ra, Hx,Hx - 1,s,p) + 1 and d^ > Aeg(ttl, M-^.h^ - 
l,s,p) + l (see [21]). ■ 

In the construction just proposed, we should choose 
to be a stronger code compared to Cx. We have given the 
construction over a nonbinary alphabet even though the case 
p = 2 might be of particular interest. 



Our next construction makes use of the cyclic finite geom- 
etry codes. Our goal will be to find a small BCH code whose 
dual is contained in a cyclic Euclidean geometry LDPC code. 
For solving this problem we need to know the cyclic structure 

of Cpp ^ (m, p., 0, s, p). Let a be a primitive element in Fpms . 

1 1 ) 

Then the roots of generator polynomial of C^q ^.(m, 0, s, p) 
are given by [11, Theorem 6], see also [13], [15]. Now, 

Z = {a^|0< max Wp.lhp'-) < (p'^ - 1 ){m - h)}, 

0<l<s 

where Wq(h.) is the q-ary weight of h. = ho + Hi q + 
••• + h-icq^^^ i.e., Wq(h.) — ^ h.|. The finite geometry 
code CgQ |,(ra, 11, 0, s,p) is actually an 1,p^) Euclidean 
geometry code. The roots of the generator polynomial of the 
dual code are given by 

Z-L={a^| min Wp. (hp'-) < h(p' - 1 )}. 

0<l<s 

In fact, the dual code is the even-like subcode of a primitive 
polynomial code of length p"^** — 1 over Fp and order m. — (o., 
whose generator polynomial, by [13, Theorem 6], has the roots 

Zp = {a^ I < min (hpM < yiiv"" - 1 )}. 

0<l<s 

Thus Z-^ = Zp U {0}. Now by [13, Theorem 11], Zp 
and therefore Z"*- contain the sequence of consecutive roots, 
a, a^, . . . , (x^°^\ where 6o = (R + ^)■p'^^ - 1 and m(p'' - 
1 ) — (m — |J.) (p'* — 1 ) = Q{v'^ — 1 ) + R. Simplifying, we see 
that R = and Q = |x giving 5o = p^'' — 1. It follows that 

4y_Jm,^,0,s,p)^ = GRMq(m,(q-l)(m-^))|Fp 
C BCn(bo]. 

Thus we have solved the problem of construction of the 
asymmetric stabilizer codes in a dual fashion to that of [10]. 
Instead of finding an LDPC code whose parity check matrix 
is contained in a given BCH code, we have found a BCH 
code whose parity check matrix is contained in a given finite 
geometry LDPC code. This gives us the following result. 

Theorem 5 (Asymmetric BCH-LDPC stabilizer codes): 
Let = C^G'j,(ra, ^,0,s,p) and 6 < 6o = P^' - 1- Let 
n = p^"^ - 1 and Cx = BCH(h] C FJ}. Then there exists an 

[[n, kx + - n, dx/dj]p 

asymmetric stabilizer code where dz > AEcfTn, \i, 1 , s, p), 

dx > 6 and kx — dim Cx, kz — dim Cz. 

Perhaps an example will be helpful at this juncture. 

Example 6: Let m = s = p = 2 and |j, = 1. Then 
Cgg' ,,(2, 1 , 0, 2, 2) is a cyclic code whose generator polyno- 
mial has roots given by 

Z = {a'^IO < max W22(2'-h) < (m- ^)(p' - 1) =3} 

0<l<2 

= {a^a^cc^a^cc^a^a^a'2}. 

As there are 4 consecutive roots and |Z| = 8, it defines a 
[15,7, > 5] code. The roots of the generator polynomial of 



the dual code are given by 

Z^ = {a^lO < min W22 (2'-h) < yiip"" - 1 ) = (2^ - 1 )} 

0<l<2 

= {a^a^a^a^a^a^a'°}. 

We see that Z^ has two consecutive roots excluding 1, 
therefore the dual code is contained in a narrowsense BCH 
code with design distance 3. Note that p*^'' — 1 =3. Thus we 
can choose Cx BCn{3) and Cz = C^y_^(2, 1 , 0, 2, 2) and 
apply Lemma|2]to construct a [[15, 3, 3/5]]2 asymmetric code. 

We can also state the above construction as in [10], that 
is given a primitive BCH code of design distance 6, find 
an LDPC code whose dual is contained in it. It must be 
pointed out that in case of asymmetric codes derived from 
LDPC codes, the asymmetry factor dx/dz is not as indicative 
of the code performance as in the case of bounded distance 
decoders. For m = p = 2, we can derive explicit relations for 
the parameters of the codes. 

Corollary 7: Let C = C^;c' c(2, 1 , 0, s, 2) and 6 = 2t + 1 < 
2** — 1 . Then there exists an 

[[22s _-i^2^'-3'- s(6 - 1),6/2'' + 1]]2 

asymmetric stabilizer code. 

Proof: The parameters of Care [2^^ — 1 ,2^*^— 3'',2^ + 1]2, 
see [15]. Since C^ is contained in a BCH code of length 
2^'^ — 1 whose design distance 6 < 2^ — 1, we can compute 
the dimension of the BCH code as 2^^ — 1 — s(6 — 1), see 
[17, Corollary 8]. By Lemma |2] the quantum code has the 
dimension 2^^^ — 3^* — s (6 — 1 ). ■ 
Example 8: For ra = p = 2 and s = 4 we can obtain a 
[255, 175, 17] LDPC code. We can choose any BCH code with 
design distance 6 < 2"* — 1 = 15 to construct an asymmetric 
code. Table I] lists possible codes. 

TABLE I 

Asymmetric BCH-LDPC stabilizer codes 



s 


5 


Code 


Asymmetry 


Rate 






[[n,k,dx/dj]2 


dz/dx 




4 


15 


[[255, 119,15/17J]2 


Ri 1 


0.467 


4 


13 


[[255, 127, 13/1 7]]2 


^ 1.25 


0.498 


4 


11 


[[255, 135, 11/1 7]]2 


1.5 


0.529 


4 


9 


[[255, 143,9/17]]2 


Ri 2 


0.561 


4 


7 


[[255, 151,7/17]]2 


^ 2.5 


0.592 


4 


5 


[[255, 159,5/17]]2 


Ri 3 


0.624 


4 


3 


[[255, 167,3/17]]2 


Ri 6 


0.655 



IV. Performance Results 

We now study the performance of the codes constructed in 
the previous section. Due to space constraints the discussion 
will be rather brief, but more details wiU be supplied in a 
forthcoming paper. We assume that the overall probability 
of error in the channel is given by p, while the individual 
probabilities of X, Y, and Z errors are Px — p/(A + 2), 
Py — P/(A+2) and Pz = pA/(A+2) respectively. The exact 
performance would require us to simulate a 4-ary channel and 
also account for the fact that some errors can be estimated 
modulo the stabilizer However, we do not account for this 



and in that sense these results provide an upper bound on the 
actual error rates. The 4-ary channel can be modeled as two 
binary symmetric channels - one modeling the bit flip channel 
and the other the phase flip channel. For exact performance, 
these two channels should be dependent, however, a good 
approximation is to model the channel as two independent 
BSCs with cross over probabilities Px + Py = 2p/(A + 2) 
and Py + Pz — p{A + 1)/(A + 2). In this case the overall 
error rate in the quantum channel is the sum of the error rates 
in the two BSCs. While this approach is going to slightly 
overestimate the error rates, nonetheless it is useful and has 
been used before [16]. Since the X-channel uses a BCH code 
and decoded using a bounded distance decoder, we can just 
compute Pg the X error rate, in closed form. The error rate in 
the Z channel, P| is obtained through simulations. The overall 
error rate is 

Pe = 1 - (1 - P^)(l - PI) = P^ + P^ - P^P| « P^ + P^. 

The LDPC code was decoded using the hard decision bit 
flipping algorithm given in [14]. The maximum number of 
iterations for decoding is set to 50. In Figure [T] we see the 
performance of [[255, 159,5/17]] as the channel asymmetry is 
varied from 1 to 100. We can clearly see the improvement as 
the channel asymmetry increases. 




Fig. 1 . Performance of [[255, 1 59, 5/1 7]] code for A = 1 , 1 0, 1 00 

The question naturally raises how do these codes compare 
with the codes proposed in [10]. Strictly speaking both con- 
structions have regimes where they can perform better than 
the other. But it appears that the algebraically constructed 
asymmetric codes have the following benefits with respect to 
the randomly constructed ones of [10]. 

> They give comparable performance and higher data rates 
with shorter lengths. 

• The benefits of classical algebraic LDPC codes are in- 
herited, giving for instance lower error floors compared 
to the random constructions. 

• The code construction is systematic. 



Our codes also offer flexibility in the rate and performance of 
the code because we can choose many possible BCH codes for 
a given finite geometry LDPC code or vice versa. The flip side 
however is that the codes given here have higher complexity 
of decoding. 
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